Periodic orbit quantization of the closed three-disk billiard as an example of a chaotic 

system with strong pruning 
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Classical chaotic systems with symbolic dynamics but strong pruning present a particular chal- 
lenge for the application of semiclassical quantization methods. In the present study we show that 
the technique of periodic orbit quantization by harmonic inversion of trace formulae, which does not 
rely on the existence of a complete symbolic dynamics or other specific properties, lends itself ideally 
to calculating semiclassical eigenvalues from periodic orbit data even in strongly pruned systems. As 
the number of periodic orbits proliferates exponentially in chaotic systems, we apply the harmonic 
inversion technique to cross-correlated periodic orbit sums, which allows us to reduce the required 
number of orbits. The power of the method is demonstrated for the closed three-disk billiard as a 
prime example of a classically chaotic bound system with strong pruning. 

PACS numbers: 05.45.-a, 03.65.Sq 



I. INTRODUCTION 

The relation between quantum spectra, on the one 
hand, and the dynamics of the corresponding classi- 
cal systems, on the other, is a problem of fundamen- 
tal importance in physics, which has attracted attention 
ever since the early days of quantum mechanics. A key 
to understanding this relation in chaotic systems was 
Gutzwiller's discovery (cf. §]) that the semiclassical 
density of states in these systems can be written as a sum 
over all (isolated) periodic orbits of the classical system. 
However, practical applications of Gutzwiller's trace for- 
mula, like other periodic orbit sums, are greatly impeded 
by the fact that the sums usually diverge in the domain 
where the physical eigenvalues are located, mainly as a 
consequence of the rapid proliferation of periodic orbits 
with growing period. Various techniques were developed 
over the years to overcome the convergence problem of 
the periodic orbit sums ^, |^, |j , but most of them depend 
on special properties of the systems, such as ergodicity 
or the existence of complete symbolic dynamics. 

As an alternative, the technique of harmonic inversion 
of semiclassical trace formulae has been introduced as 
a very general method for extracting eigenvalues from 
periodic orbit sums for both chaotic and regular sys- 
tems [|[ 1^, |8[ ||, |l^ . Harmonic inversion circumvents 
the convergence problems, and allows one to obtain the 
semiclassical eigenvalues from a relatively small number 
of periodic orbits. Moreover, the method does not re- 
quire any special properties of the systems, and therefore 
is applicable also to the particularly challenging case of 
systems with strong pruning, i.e., systems with a highly 
incomplete symbolic coding of the periodic orbits. In 
this paper we demonstrate the power of the harmonic 
inversion technique in the semiclassical quantization of 
the closed three-disk billiard system, for which all other 
semiclassical quantization methods to date have failed in 
determining more than the few lowest eigenvalues, chiefly 
because of the extremely rapid proliferation of periodic 
orbits with increasing action and the strong pruning typ- 



ical of this system. From a technical point of view our 
study will show that harmonic inversion is numerically 
stable even when handling huge (on the order of several 
million) periodic orbit sets. 

The three-disk billiard system consists of three equally 
spaced hard disks of unit radius. The existence or nonex- 
istence of periodic orbits and the behavior of the periodic 
orbit parameters in the three-disk system sensitively de- 
pend on the distance d between the centers of the disks. 
Large disk separations, especially d = 6, have served as 
a test case for periodic orbit quantization of chaotic sys- 
tems in many investigations in recent years. In particu- 
lar, the system has been employed as prototype example 
for the usefulness of periodic orbit quantization by cycle 
expansion techniques |3[ |2, 13, 1^. For large disk sepa- 
rations, the assumptions of the cycle expansion, namely 
that the contributions from long orbits are shadowed by 
those of short orbits, are almost exactly fulfilled. This is 
no longer true for small disk separations. As the disks 
approach each other, the convergence of the cycle expan- 
sion becomes slower and slower, until it finally breaks 
down. 

Periodic orbit quantization by harmonic inversion does 
not depend on specific properties, such as shadowing of 
orbits etc., and therefore is expected to work well also 
for small disk separations. In fact, the harmonic inver- 
sion technique has successfully been applied to the "stan- 
dard" literature disk separation d = 6 as well as to the 
small distance d = 2.5, where the convergence of the cycle 
expansion is already slow 0, 1). In this paper we con- 
centrate on the limiting case of mutually touching disks, 
d — 2, where the system turns into a bound system, and 
the requirements of the cycle expansion and other semi- 
classical methods are no longer fulfilled at all. 

For the closed three-disk system, the lowest exact 
quantum eigenvalues were first calculated by Tanner et 
al. and by Scherer A remarkable step towards 
the semiclassical quantization of the closed three-disk 
system was achieved by combining the conventional cy- 
cle expansion with a functional equation Based on 
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periodic orbit corrections to the mean density of states, 
approximations to the very lowest eigenvalues could be 
obtained from a small set of periodic orbits up to cy- 
cle length ^ = 3, where the symbolic dynamics is still 
complete (including formally a zero length orbit). For 
the very lowest eigenvalues, the results were found to be 
in good agreement with the exact quantum values. How- 
ever, the method failed for higher eigenvalues, where only 
the mean density of states could be reproduced. Since the 
problems with the method arose from the strong prun- 
ing of orbits in this system, they are of a fundamental 
nature, and could not be overcome by including more or- 
bits. No other semiclassical method has been successful 
in obtaining higher eigenvalues in the closed three-disk 
billiard system yet. In this paper, we will demonstrate 
that harmonic inversion passes even this acid test of semi- 
classical quantization. 

In Sec. II we discuss the construction of the periodic 
orbit signal for the closed three-disk system, which is a 
nontrivial task because of, firstly, the extremely rapid 
proliferation of periodic orbits and, secondly, the strong 
pruning of the symbolic dynamics. Once the signal has 
been constructed, we refer in Sec. Ill to the harmonic 



inversion technique for the extraction of the semiclassi- 
cal eigenvalues. The numerical recipe of the harmonic 
inversion method has already been established in the lit- 
erature 1^, 0, H, H, |l^, |ri|, however, to make the ar- 
ticle se lf-con sistent we recapitulate the important steps 
in Sec. Ill A| . The s emiclassical results are discussed in 
Sees. |III B| and pT^ . 



II. PERIODIC ORBITS OF THE CLOSED 
THREE-DISK SYSTEM 

We start by summarizing the basic properties of the 
three-disk system relevant for periodic orbit quantiza- 
tion. Furthermore, the pruning of orbits and the distri- 
bution of the orbit parameters are investigated, and the 
resulting strategies for the numerical search for periodic 
orbits are discussed. 



A. Symbolic code and symmetry reduction 

The periodic orbits of the three-disk system can be 
labeled by a ternary symbolic code, which is complete 
for sufhciently large disk separations. The classical dy- 
namics of the three-disk system has been studied, among 
others, by Gaspard and Rice If the disks are labeled 
by the numbers 1, 2, 3, each periodic orbit is charac- 
terized by a sequence of these numbers, indicating the 
disks the particle collides with during one period of the 
orbit. For sufficiently separated disks, there is a one-to- 
one correspondence between the symbolic code and the 
periodic orbits of the system: For every sequence there 
exists one unique periodic orbit (with the restriction that 
consecutive repetitions of the same symbol are forbidden 
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FIG. 1: Open and closed three-disk system. The fundamen- 
tal domain consists of a one-sixth slice of the full system. 



and circular shifts of a sequence describe the same or- 
bit). As was first pointed out by Hansen at disk 
separation d = 2.04821419 "pruning" sets in, i.e., part 
of the periodic orbits become unphysical as they begin 
to run through one of the disks, and thus the symbolic 
code is no longer complete. The number of pruned orbits 
rapidly becomes larger and larger if the disks continue 
to approach each other. In the limiting case of touching 
disks, d = 2, the system exhibits strong pruning. Exam- 
ples of pruned orbits will be discussed in Section [H C| . 

In our calculations, we make use of the symmetry 
reduction introduced by Cvitanovic and Eckhardt in 
Ref. 1^: The three-disk system is invariant under the 
symmetry operations of the group C^y, i.e., reflections at 
three symmetry lines and rotations by 27r/3 and 47r/3. 
The periodic orbits fall into three classes of distinct sym- 
metry: orbits invariant under reflections at one of the 
symmetry lines (multiplicity 3), orbits invariant under 
rotations by 27r/3 and 47r/3 (multiphcity 2), and orbits 
with no symmetry (multiplicity 6). The quantum states 
are grouped in the three irreducible subspaces Ai, A2, 
and E, where the states of the Ai {A2) subspace are sym- 
metric (antisymmetric) under reflection at the symmetry 
lines, respectively, and the states of the E subspace are 
invariant under rotations by 27r/3 and 47r/3. 

Following Refs. [|,|l|, one can map the system onto a 
fundamental domain, which consists of a one-sixth slice 
of the full system, with the symmetry axes acting as 
straight mirror walls (see Fig. |l|). The periodic orbits 
of the full system can be described completely in terms 
of the periodic orbits in the fundamental domain. The 
symmetry reduced periodic orbits are labeled by a binary 
code, where the symbol '0' represents backscattering (or 
change between clockwise and anti-clockwise scattering) 
and the symbol '1' stands for scattering to the third disk 
in the original full domain picture. For example, the 
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shortest full domain orbit 12 maps onto the orbit, and 
the 123 orbit maps onto the 1 orbit. 



B. Semiclassical density of states 

The quantum resonances of the three different sub- 
spaces Ai , A2 and E can be obtained separately from the 
periodic orbits in the fundamental domain by introducing 
appropriate weight factors for the orbits in Gutzwiller's 
trace formula, as was shown by Cvitanovic and Eckhardt 
pof . We concentrate on the Ai subspace, for which each 
orbit has a weight factor equal to 1. 

As for all billiard systems, the shape of the perio dic 
orbits is independent of the wave number k = \^2mE/h, 
and the action scales as 



S/h = ks, 



(1) 



where the scaled action s is equal to the physical length of 
the orbit. We consider the density of states as a function 
of the wave number 



p{k) = - llm g{k) , 

TT 



(2) 



with a scaled response function g{k). Applying the 
Gutzwiller trace formula to the closed three-disk system 
yields for the Ai subspace 



7(fc) = <?o(fc)-iE 
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r|det(Mpo-l)|i/2 
r|(Apo-l)(^-l)|V2 



(3) 



where Mpo and are the monodromy matrix and the 
Maslov index of the orbit, respectively. Is is the symbol 
length, Spo is the scaled action, and Apo denotes the ex- 
panding stability eigenvalue of the orbit (i.e., the eigen- 
value with an absolute value larger than one). The sum 
runs over all symmetry reduced periodic orbits including 
multiple traversals. Here, r denotes the repetition num- 
ber with respect to the corresponding primitive orbit. 

In practice, only the primitive periodic orbits have to 
be determined. The parameters of the r'^ repetition of 
the primitive orbit (here characterized by the index 0) 
are then given by Is — rl so, s = rso and A = Aq. 



C. Numerical search for periodic orbits 

For extracting the quantum resonances of the system 
from Eq. (Q) by harmonic inversion, all periodic orbits up 
to a maximum scaled action have to be included. The pa- 
rameters of the periodic orbits - scaled action and stabil- 
ity eigenvalues - have to be determined numerically. We 
calculate the primitive periodic orbits using the symbolic 
code as input. For simplicity, the calculations are carried 





FIG. 2: An example of the first class of pruned orbits (see 
text), plotted in full domain representation. The symmetry 
reduced code of the orbit is 000011. The orbit is shown at 
two different disk separations d, indicated at the top of each 
diagram. At d = 2.04, the orbit is still physical, at d = 2.0 it 
has become unphysical as is penetrates the disks. 



out in the full domain, and the results are then trans- 
lated back into the symmetry reduced system. The disks 
are "connected" according to the code, starting with ar- 
bitrary reflection points on the disks as initial condition. 
The reflection points are then varied in such a way that 
the total length (i.e., the action) of the orbit reaches a 
minimum. All orbit parameters can then be calculated 
from the reflection points. The scaled action is given by 
the length of the (symmetry reduced) orbit, and the sta- 
bility eigenvalue A can be determined by an algorithm 
proposed by Bogomolny |2^. In addition to action and 
stability, one can also determine averages of different clas- 
sical quantities (distance from the center of the system, 
angular momentum, etc.), which are needed for the cross- 
correlation technique (see Section [II). 



For disk separations smaller than the pruning limit 
d — 2.04821419, it has to be checked whether or not the 
orbits are physical, i.e., whether or not they stay com- 
pletely outside the disks. In our numerical calculations, 
we found two different classes of pruned orbits. An ex- 
ample of the flrst class is given in Figure ^: As the disks 
approach each other, a section of the orbit connecting 
two of the disks gets inside the third disk. The pruned 
orbit still corresponds to a unique minimum of the total 
length when the reflection points on the disks are varied 
for given symbolic code. This is not the case for the sec- 
ond class, an example of which is shown in Figure |^: In 
this case, as the disks come closer, the reflection angle at 
one specific reflection point approaches n. As the value tt 
is reached, the minimum of the action with respect to this 
specific reflection point splits into two equal minima and 
one local maximum. Accordingly, the orbit splits into 
three unphysical orbits with the same symbolic code. 

In one case, the orbit is now reflected on the inside of 
the disk (see Fig. ^), which corresponds to a reflection 
angle larger than tt (i.e., the orbit has penetrated into 
the disk in the vicinity of the reflection point). This or- 
bit does not correspond to a minimum of the total length 
but to a saddle point. In the other two cases. Figs. ||d 
and the reflection point moves along the disk in such 
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FIG. 3: An example of the second class of pruned orbits (see 
text), plotted in full domain representation. The symmetry 
reduced code of the orbit is 0000001. The orbit is still physical 
at disk separations d = 2.04 and d — 2.02. At d = 2.0 it has 
split into three unphysical orbits. The squares indicate the 
"reflection points". Note in (c) that the reflection angle at 
the three central reflection points is larger than tt, i.e., the 
orbit in fact penetrates into the disks in the vicinity of these 
points. 

a way that the reflection law is no longer fulfilled but the 
orbit just passes the "reflection point" in a straight line 
and penetrates into the disk. These two cases correspond 
to two equally deep minima of the total length. The two 
orbits have the same shape but only differ as to which 
penetration point is considered as "reflection point" , in- 
dicated by the squares in Figs. ||d and In fact, the 
shape of the orbit is equal to that of a pruned orbit with 
a (symmetry reduced) symbol length that is shorter by 
1 and which belongs to the first class of pruned orbits 
as defined above. (In this orbit the reflection point in 
question is simply missing, i.e., none of the penetration 
points is considered as reflection point.) 

Besides the existence or nonexistence of periodic or- 
bits, the distance between the disks also strongly influ- 
ences the distribution of the periodic orbit parameters. 
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FIG. 4: Distribution of periodic orbit parameters of the 
three-disk system (a) with disk separation d = 6, and (b) in 
the case of touching disks, d = 2. The quantities plotted are 
the absolute value of the larger stability eigenvalue |A| versus 
the physical length s of the orbits. For d — 2, for clarity a few 
channels of orbits are marked by arrows. Channels (A) and 

(B) break off because of pruning. The two channels marked 

(C) contain infinitely many orbits, but have been cut off by 
restricting the maximum number of successive zeros in the 
symbolic code (see text). 

This is illustrated in Figure |^, which compares the or- 
bit parameters of the shortest primitive periodic orbits 
found for the open three-disk system with disk separa- 
tion d — 6 and for the case of touching disks, d = 2. For 
d — 6, the action and stability of the orbits are mainly 
determined by the symbol length, the orbit parameters 
are nicely bunched in easily distinguishable intervals, and 
the number of orbits up to a given action is relatively 
small. However, the picture changes completely when the 
disks approach each other: As all orbits become shorter, 
the total number of orbits up to a given action increases 
rapidly. The parameters of the orbits are no longer sim- 
ply determined by the symbol length, and their overall 
behavior becomes much more complicated, as is evident 
from Figure ^. 

For d = 2, the orbits can be grouped in "channels" with 
the same "tail" (end figures) but growing number of lead- 
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FIG. 5: Examples of orbits of the 0"101 channel, (a) - (c), 
and the 0"111 channel, (d) - (f) (dubbed (C) in Fig. gb), in 
full domain representation. With growing number of leading 
'O's, the orbits run deeper and deeper into the corners between 
the disks. The series converges towards the limiting orbit that 
starts exactly in the point where the disks touch and is then 
directly reflected back from the opposite disk. 



ing 'O's in the code. (A sequence of n leading 'O's in the 
code will be denoted by 0" in the following). These or- 
bits have the same basic shape but run deeper and deeper 
into the corner formed by two touching disks, bouncing 
back and forth between the two disks (see Fig. ||). In 
each channel, the action of the orbits grows very slowly 
with increasing symbol length, while the expanding sta- 
bility eigenvalue increases relatively fast (in fact expo- 
nentially) with every additional collision. Since adding 
a leading '0' to the symbolic code does not change the 
action considerably, there is a huge number of orbits with 
very long symbol lengths but relatively small, and com- 
parable, actions. In fact, most channels break off because 
of pruning. E.g., the 0"! series, labeled (A) in Fig. 
breaks off after six orbits, and the 0"11 series, marked 
(B), already breaks off after three orbits. Our analysis 
strongly suggests that there exist only two infinite chan- 



nels: The 0"101 and the 0"111 series, dubbed (C) in 
Fig. ||b (these channels are only shown up to 20 leading 
'O's in the code). The shapes and actions of these orbits 
converge very slowly towards the shape and action of the 
limiting orbit which starts exactly in the touching point 
of two disks and is directly reflected back by the oppo- 
site disk (cf. Fig. ||). Although all other channels break 
off because of pruning, the number of orbits in a chan- 
nel up to a given action often becomes very large (or, in 
the two cases mentioned, even infinite). Note that the 
orbit, bouncing back and forth between two disks, does 
not exist in the closed three-disk system. 

While for large disk separations finding all orbits up 
to a given action simply corresponds to calculating all 
orbits up to a certain maximum symbol length, the com- 
plicated distribution of periodic orbit parameters in the 
closed three-disk system renders the search for the rel- 
evant periodic orbits up to a given action a nontrivial 
task. Due to the extremely rapid proliferation of peri- 
odic orbits with increasing length s, we were not able 
to determine all orbits up to a given length Smax, but we 
had to introduce restrictions for the orbits to be included. 
As the amplitude in the Gutzwiller formula depends on 
the inverse of the stability parameter, a reasonable cut- 
off criterion is the size of the stability of the orbit. On 
the other hand, since in each channel the absolute value 
of the expanding stability eigenvalue grows exponentially 
fast, the main contributions of each channel to the pe- 
riodic orbit sum come from the orbits with small sym- 
bol lengths. Furthermore, the Maslov indices of adjacent 
orbits in each channel differ by two, and therefore the 
periodic orbit amplitudes alternate in sign, which leads 
to approximate cancellations of terms in the semiclassi- 
cal trace formula. As a second criterion, we therefore 
restricted the number of consecutive symbols '0' in the 
symbolic code, thus selecting the most relevant orbits of 
each channel. 

The orbits were calculated channel by channel, start- 
ing from the shortest orbit in the channel and adding 
more and more leading zeros to the code. The calcula- 
tion of the channel was broken off if one of the following 
conditions was fulfilled: 

• The action exceeds the maximum action Smax- 

• The absolute value of the stability eigenvalue ex- 
ceeds the given value Amax- 

• The maximum number of consecutive symbols '0' 
is reached. 

• The orbits become pruned. 

For the set of periodic orbits shown in Figure we 
searched for orbits with physical length s < Smax = 5.0. 
The maximum for the absolute value of the stability 
eigenvalue was chosen to be Amax = 10^. The maximum 
number of successive zeros in the code was restricted to 
20 for orbits with length s < 3.1 and to 12 for orbits with 
3.K s < 5.0. 
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III. PERIODIC ORBIT QUANTIZATION OF 
THE CLOSED THREE-DISK SYSTEM BY 
HARMONIC INVERSION 

The fundamental problem, and challenge, now is to 
extract semiclassical eigenvalues in a numerically stable 
way from the huge periodic orbit set presented in Fig. 
To this end, we resort to harmonic inversion of periodic 
orbit signals 0, |[ ^, |l^ [T^. Here, we review the basic 
ideas and present an extension of the method introduced 
in pH l to harmonic inversion of cross-correlated periodic 
orbit sums. 



A. The method 

The starting point is to introduce a weighted density 
of states in terms of k 



Qaa' {k) 



1 



with 



(4) 



(5) 



where km is the eigenvalue of the wave number of eigen- 
state |m) and 



{m\Aa\m) 



(6) 



are the diagonal matrix elements of a chosen set of D 
linearly independent operators Aa, a — 1, 2, . . . , D. The 
Fourier transform of (||) yields a D x D cross-correlated 
signal 
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+ 00 
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(7) 



A semiclassical approximation to the weighted quantum 
response function is given by the semiclassical re- 
sponse function weighted with the classical averages 



5po Jo 



^a(q(s),p(s))ds 



(8) 



with Aa (q, p) the Wigner transform of the operator A^ 
pl| , . By Fourier transformation we obtain the semi- 
classical approximation to the quantum signal (0), which 
reads 



po 



^Q,po O!'a',po -^po 6 2 f^p° 

r|det(Mpo-l)|i/2 



5{s- Spo) 



(9) 



where r is the repetition number counting the traversals 
of the primitive orbit, and Mpo and /Xpo are the mon- 
odromy matrix and Maslov index of the orbit, respec- 
tively. Semiclassical approximations to the eigenvalues 
km, and eventually also to the diagonal matrix elements 
{m\Aa\m), are obtained by adjusting the semiclassical 
cross-correlated periodic orbit signal (M to the functional 
form of the quantum signal (|^) . This can be done with 
the filter-diagonalization method for the signal process- 
ing of cross-correlation functions ||2^, U|. Here, we resort 
to an alternative method introduced in Ref. [ pl] |, which 
will now be extended for the harmonic inversion of cross- 
correlated periodic orbit sums. 

The special form of the periodic orbit signal (|9|) as a 
sum of 5 functions allows for an even simpler procedure 
than filter-diagonalization, viz. analytical filtering. In the 
following we will apply a rectangular filter, i.e., f{k) = 1 
for frequencies k G [ko — Ak,ko + Ak], and f(k) — 
outside the window. Applying this filter to the semiclas- 
sical signal C^^, (s) in (|^) we obtain the band-limited (bl) 
cross-correlated periodic orbit signal, 



Jko-Ak 
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(10) 



5poy 



The introduction of ko into the arguments of the expo- 
nential functions in causes a shift of frequencies by 
— fco in the frequency domain. Note that C'^^^'(s) is a 
smooth function and can be easily evaluated on an ar- 
bitrary grid of points Sm < Smax provided the periodic 
orbit data are known for the set of orbits with classical 

action Spo < Smax- 

Applying now the same filter as used for the semiclas- 
sical periodic orbit signal to the quantum one, we obtain 
the band-limited quantum signal 

f fco+Afc 



feo-Afc 



M 

E 

rn— 1 



-i{kjn—kQ)s 



fcol < Afc . (11) 



In contrast to the signal C^^,{s) in Eq. (0), the band- 
limited quantum signal consists of a finite number of 
frequencies km, = 1,...,M, where in practical ap- 
plications M can be of the order of ~ (50-200) for an 
appropriately chosen frequency window, Afc. The prob- 
lem of adjusting the band-limited semiclassical signal in 
Eq. (^^ to its quantum mechanical analogue in Eq. (|TT 
can now be written as a set of nonlinear equations 



n = 0,l,...,2N-l , 



po 



M 



, (12) 
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for the unknown variables, viz. the shifted frequencies, 
k'^ = km — kQ, and the parameters bam, m = 1, . . . , M. 
Note that M is related to the signal length 2N and the 
dimension D of the cross-correlation matrix by M = ND. 
The signal (^2|) can be evaluated on an equidistant grid, 
s = nr, with relatively large step size t = tt/Aw, and 
the discrete signal points Caa',n = C^^^{nT) are called 
the "band-limited" cross-correlated periodic orbit signal. 

We now wish to solve the nonlinear system, Eq. (12), 
which can be written as 



M 

E 

m— 1 



(13) 



where Zm = exp (— ifc^r) and bam are, generally com- 
plex, variational parameters. We assume that the num- 
ber of frequencies in the signal is relatively small {M ^ 50 
to 200). Although the system of nonlinear equations 
is, in general, still ill-conditioned, frequency filtering re- 
duces the number of signal points, and hence the num- 
ber of equations. In Ref. |ll| three different methods, 
viz. linear predictor (LP), Pade approximant (PA), and 
signal diagonalization (SD) have been employed to solve 
Eq. ( p^ for a one-dimensional signal, i.e., D = 1. The 
SD method can be generalized in a straightforward man- 
ner to the harmonic inversion of cross-correlated periodic 
orbit sums. The problem of solving the nonlinear set of 
equations ( |l3|) can be recast in the form of the general- 
ized eigenvalue problem [Eil p5l p6[ 
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(14) 



where the elements of the M x M operator matrix U and 
overlap matrix S depend trivially upon the signal points: 



Caa'.i+j] i, j = 0, . . . ,7V- 1 . (15) 



The matrices U and S in Eq. ( p^ are complex symmet- 
ric (i.e., non-Hermitian) , and the eigenvectors Bm are 
orthogonal with respect to the overlap matrix S, 



{Bm\S\Bm') 



N () , 



(16) 



where the brackets define a complex symmetric inner 
product {a\b) = (6|a), i.e., no complex conjugation of ei- 
ther a or b. The overlap matrix S is not usually positive 
definite and therefore the A'^'s are, in general complex, 
normalization parameters. An eigenvector Bm cannot be 
normalized for Nm = 0. The parameters bam in Eq- ^ 
are obtained from the eigenvectors Bm via 



(17) 



The parameters Zm in Eq. (|13|) are given as the eigen- 
values of the generalized eigenvalue problem ( |l4|), and 
are simply related to the frequencies k'm in Eq. (|12D via 
Zm = exp(-ifc^r). 




The resolution of the results depends on the signal 
length. For a one-dimensional signal, the method re- 
quires a signal length of Smax ~ 47rp(fc), with 3(k) the 
mean density of states, to resolve the frequencies g . This 
means that all periodic orbits up to the scaled action 
•Smax have to be included. The advantage of using the 
cross-correlation approach is based on the insight that 
the total amount of independent information contained 
in the D x D signal is D{D + 1) multiplied by the length 
of the signal, while the total number of unknowns (here 
bam and km) is {D + 1) times the total number of poles 
km- Therefore the informational content of the D x D 
signal per unknown parameter is increased (compared to 
the one-dimensional signal) by roughly a factor of D, and 
the cross-correlation approach should lead to a significant 
improvement of the resolution. The power of this method 
for periodic orbit quantization has already been proven 
for the example of an intcgrable system We will 

now demonstrate how the method works in the case of 
three touching disks, as an example of a chaotic system 
with extremely rapid proliferation of orbits and strong 
pruning. 



B. Semiclassical eigenvalues of the closed 
three-disk system by harmonic inversion of a 
one-dimensional signal 

We first investigate a one-dimensional signal (A^ =1) 
by simply choosing Ai = 1, i.e., the unity operator. The 
signal is constructed from the set of periodic orbits shown 
in Fig. which contains about 5 million primitive or- 
bits. 

We analyzed the signal up to length s„iax ==_4^ using 

The 



the decimated signal diagonalization method 1 1 1 



frequencies obtained by harmonic inversion of the signal 
are the semiclassical approximations to the eigenvalues 
kn- Fig. @ shows the results from harmonic inversion 
(solid lines), compared with the exact quantum results 
(dashed lines), both presented in terms of the energy 
E = (Rek)'^/2. The exact quantum eigenvalues were 
taken from Scherer and Wirzba |2^. The harmonic 
inversion results clearly reproduce the quantum eigenval- 
ues up to £^ « 4500 with only small deviations. However, 
the signal length was not sufficient to resolve the eigen- 
values in the region E > 4500, where the density of states 
with respect to the wave number k (which grows ^ \/E) 
becomes too large. 

Apart from the values shown in Fig. |^, the harmonic 
inversion of the signal also yielded a number of uncon- 
verged frequencies. Converged frequencies were identi- 
fied by having an imaginary part close to zero (since the 
eigenvalues of the wave number must be real for bound 
systems) and an amplitude close to the theoretical value 
mfe = 1. The whole set of frequencies obtained is shown 
in Figure |^ together with the real part of their ampli- 
tudes. The solid vertical lines mark the positions of the 
exact quantum eigenvalues. The semiclassical values in- 
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FIG. 6: Energy eigenvalues of the closed three-disk system 
(j4i subspace). Dashed lines: exact quantum eigenvalues; 
solid lines: results E = (Re fc)^/2 from harmonic inversion 
of a signal of length Smax = 4.9. 
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FIG. 7: Amplitudes of the whole set of values E = (Re fe)^/2 
(all circles) obtained for the energy eigenvalues of the closed 
three-disk system by harmonic inversion (including also un- 
convereed values) . The filled circles mark the values presented 
in Fig. |6[ The solid vertical lines indicate the positions of the 
exact quantum eigenvalues. 



eluded in Fig. |^ are represented by filled circles in Fig. 
The amplitudes are not as well converged as the frequen- 
cies (i.e. the eigenvalues) but partly show larger devia- 
tions from the theoretical value ruk = 1. [This can be 
explained by the observation that when the decimated 
signal diagonalization method is used for harmonic in- 
version the frequencies usually converge faster than the 
amplitudes.] However, up to « 4500, the eigenvalues 
can be clearly identified. 

The small deviations of the eigenvalues obtained by 
harmonic inversion from the exact quantum values may 
have different reasons. Especially for the lowest eigenval- 
ues, the major part is probably due to the semiclassical 
error. On the other hand, as explained above, we did not 




7000 



FIG. 8: Harmonic inversion results E = {Rekf /2 for the 
energy eigenvalues of the closed three-disk system as a func- 
tion of the signal length Smax. Only results with |Imfc| < 1.0 
and amplitudes Rem^ > 0.2 are included. The filled circles 
mark values with Rem^ > 0.5. (Solid vertical lines: positions 
of the exact quantum energies.) 



include all orbits in the signal but left out a large num- 
ber of relatively unstable orbits. Although each of these 
orbits only gives a negligible small contribution to the 
periodic orbit sum, the number of excluded orbits may 
be so large that their contributions sum up in a way so as 
to have a visible effect on the density of states. Finally, 
even in the low-lying part of the spectrum, deviations 
may arise from the fact that the signal length was very 
short. The influence of the missing orbits and the rel- 
atively short signal length are reflected in the relatively 
poor convergence of the amplitudes. 

In order to test the stability of the results with respect 
to the signal parameters, and to obtain an estimate of 
random errors, we performed the same calculation with 
various sets of different signal parameters. Figure ^ shows 
the harmonic inversion results for the eigenvalues as a 
function of the signal length Smax- From all frequencies 
obtained the values with jlm/cj < 1.0 and with ampli- 
tudes Remfe > 0.2 were singled out, respectively. The 
filled circles mark amplitudes with Remfc > 0.5. It 
can clearly be seen how the maximum eigenvalue up to 
which the spectrum can be resolved depends on the sig- 
nal length. Fig. ^ shows that the results for the low-lying 
eigenvalues are quite stable with respect to variation of 
the signal length. 

Compared with the results obtained by Tanner et 
al. |]l5|] using the extended cycle expansion method, the 
analysis of a single signal by harmonic inversion has 
yielded about the same number of eigenvalues for the 
closed three-disk system. For the resolution of higher 
eigenvalues, the signal length was not sufficient. Extend- 
ing the signal to significantly larger scaled actions is in 
practice not possible because of the extremely rapid in- 
crease of the number of orbits. It should be stressed. 
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FIG. 9: Energy eigenvalues of the closed three-disk sys- 
tem [Ai subspace). Dashed lines: exact quantum eigenvalues; 
solid lines: results E = (Refc)^/2 from harmonic inversion of 
a 3 X 3 cross-correlated signal with length Smax = 4.8. The op- 
erators chosen for the construction of the signal are 1 (unity), 
r* and L*, where r and L are the distance from the center of 
the system and the absolute value of the angular momentum, 
respectively. 



however, that, in contrast to the cycle expansion method, 
which runs into severe problems because of the pruning 
of orbits in this system, it is at least in principle possible 
to improve the harmonic inversion results by including 
more orbits. On the other hand, as discussed before, the 
harmonic inversion method offers the possibility to sig- 
nificantly reduce the signal length required for the reso- 
lution of eigenvalues by the construction and analysis of 
cross-correlated periodic orbit sums. 



C. Improvement of the resolution by harmonic 
inversion of cross-correlated periodic orbit sums 

To increase the accuracy of the semiclassical eigenval- 
ues, we now apply the cross-correlation technique to the 
closed three-disk system. Instead of the one-dimensional 
signal investigated in the previous section, we choose a set 
of three operators and construct a 3 x 3 cross-correlated 
signal according to Eq. (^). As input we again use the 
set of periodic orbits presented in Fig. ^ together with the 
averages of different classical quantities over the orbits. 

Fig. ^ shows the results from the harmonic inversion 
of a 3 X 3 signal of length Smax — 4.8 (soHd lines), which 
was constructed using the operators 1 (unity), and L**, 
where r and L denote the distance from the center of the 
system and the absolute value of the angular momen- 
tum, respectively. The eigenvalues are again presented 
in terms of the energy E — {KekY /2. For comparison, 
the dashed lines indicate the positions of the exact quan- 
tum eigenvalues. The converged eigenvalues for the wave 
number k have been singled out from the whole set of fre- 
quencies obtained by the condition that they should have 
an imaginary part close to zero and an amplitude close 
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FIG. 10: Closed three-disk system: The whole set of fre- 
quencies (all circles) resulting from the analysis of the 3x3 
signal (see Fig. ^). (a) Real part of the amplitudes versus 
energy E = (Refc)^/2. The solid vertical lines mark the po- 
sitions of the exact quantum eigenvalues, (b) Imaginary part 
of the frequencies versus real part of the amplitudes. In both 
diagrams, the filled circles designate the values included in 
Fig. |. 



to the theoretical value mk = 1. The amplitudes and the 
imaginary parts of all frequencies obtained are shown in 
Figure |l^ by circles. In particular, the filled circles de- 
note the frequencies included in Fig. |[ In Fig. the 
positions of the exact quantum eigenvalues are given by 
the solid vertical lines. 

With the cross-correlated signal, one can now clearly 
identify eigenvalues up to ~ 6500. In addition, the 
convergence of the lowest eigenvalues is improved as com- 
pared t o the results from the single signal obtained in 
Section IIIB , as can be seen from the amplitudes. It is 
not possible to determine to what extent the accuracy 
of the semiclassical eigenvalues has improved since the 
results can only be compared with the exact quantum 
eigenvalues and the size of the semiclassical error is un- 
known. 
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FIG. 11: Harmonic inversion results E = (Refc)^/2 for the 
energy eigenvalues of the closed three-disk system as a func- 
tion of the signal length Smax, calculated (a) from a 3 x 3 
cross-correlated signal including the operators 1 (unity), r** 
and I/'* and (b) from a 4 x 4 signal including the operators 1, 
, and (r: distance from the center of the system, L: 
absolute values of the angular momentum) . Only results with 
|Imfc| < 1.0 and amplitudes Remfc > 0.5 are included. The 
filled circles mark values with |Imfc| < 0.5. The solid vertical 
lines indicate the positions of the exact quantum eigenvalues. 



Again, we have performed the same calculation for var- 
ious sets of parameters in the harmonic inversion scheme. 
Figure |ll| shows the results from a 3 x 3 and from a 
4x4 signal as functions of the signal length. The 3x3 
signal was constructed from the operators 1 (unity), 
and L^, and the 4x4 signal contained the operators 1, 
r^, r** and L^. In both diagrams, only the frequencies 
with an imaginary part |Imfc| < 1.0 and an amplitude 
Remfc > 0.5 were included. The filled circles mark values 
with |Imfc| < 0.5. The positions of the exact quantum 
eigenvalues are again marked by the solid vertical lines. 

As with the one-dimensional signal analyzed in the pre- 
vious section, it can be observed that the results for the 
lowest eigenvalues are very stable with respect to the vari- 



ation of the signal length as long as the signal is not too 
short. However, the very lowest frequencies show a ten- 
dency to split into two or more, which was a frequent ob- 
servation in our calculations when the signal length and 
the matrix dimensions were chosen too large. The results 
for the higher eigenvalues depend more sensitively on the 
signal length as the density of states approaches the lim- 
iting resolution that can be achieved with the signal. In 
this region, the 4x4 signal shows a better resolution than 
the 3x3 signal. 

With the 3x3 signal and the two largest values of 
the signal length considered, the lowest eigenvalue was 
not obtained. The reason for this probably lies in in- 
accuracies at the end of the signal, as the signal length 
approaches the value Smax = 5.0: In the decimated sig- 
nal diagonalization scheme, which we used for the har- 
monic inversion of the signal, the 5 functions in the 
original signal (^) are replaced with sine functions jl^ 
[sine X = (sin2;)/a;]. Therefore, every periodic orbit con- 
tributes to the signal not only exactly at its scaled ac- 
tion s — Spo, but also in a range As around this point. 
If the signal length is close to Smax = 5.0, also orbits 
with actions slightly larger than this limit may give non- 
negligible contributions to the very end of the signal. 
Since we did only include orbits with s < 5.0, these con- 
tributions are missing in our signal. 

Some eigenvalues between E = 3500 and E ~ 5000 
seem to be particularly hard to obtain from the set of 
orbits used. Possibly, these eigenvalues are related to the 
orbits running very deep into the corners between the 
disks, which were not included in the signal. We conjec- 
ture that the missing orbits as well as the short signal 
length are again responsible for the deviations especially 
of the results for higher eigenvalues from the exact quan- 
tum values. 



IV. CONCLUSION 

We have successfully applied the technique of harmonic 
inversion to the semiclassical quantization of a system 
with an extremely rapid proliferation of periodic orbits 
and strong pruning. With the harmonic inversion of 
cross-correlated periodic orbit sums we were able to cal- 
culate accurate semiclassical energy eigenvalues of the 
closed three-disk system up to the region E « 6500 (a 
total of 29 levels). This is the first time that such high- 
lying eigenvalues have been correctly calculated by peri- 
odic orbit theory in this system. Our results demonstrate 
that the harmonic inversion of cross-correlated periodic 
orbit sums is indeed a powerful method of periodic or- 
bit quantization and allows the stable handling even of 
huge periodic orbit sets. In contrast to other methods, 
which depend on the existence of a complete symbolic 
code, the harmonic inversion technique is not affected 
by the strong pruning of orbits. In the closed three-disk 
system, the only restriction for the practical application 
of the harmonic inversion method lies in the extremely 
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large number of periodic orbits, which is a special feature 
of this system. It should also be noted that the semiclas- 
sical eigenvalues were obtained without resorting to the 
mean density of states, which, including all necessary cor- 
rection terms, is nontrivial to calculate even in the case 
of the closed three-disk system . The present results 
may stimulate future work on other challenging systems 
without special structural information on the quantum 
or classical level. One example with strong pruning of 
the symbolic dynamics is the three-dimensional general- 
ization of the closed three-disk billiard, i.e., a system con- 
sisting of four touching spheres at the corners of a regular 
tetrahedron. By contrast to the closed three-disk system 



this is an open system where no functional equation can 
be applied. Experiments on chaotic light-scattering from 
the four spheres have recently attracted much attention 
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